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Using previously developed exciton scattering model accounting for the interband, i.e., exciton- 
biexciton, Coulomb interactions in semiconductor nanocrystals (NCs), we derive a closed set of 
equations for 2D double- quantum coherence signal. The signal depends on the Liouville space path- 
ways which include both the interband scattering processes and the inter- and intraband optical 
transitions. These processes correspond to the formation of different cross-peaks in the 2D spectra. 
We further report on our numerical calculations of the 2D signal using reduced level scheme param- 
eterized for PbSe NCs. Two different NC excitation regimes considered and unique spectroscopic 
features associated with the interband Coulomb interactions are identified. 



I. INTRODUCTION 

Coulomb interactions between carriers in semiconduc- 
tor materials including nanocrystals (NCs) play central 
role in photoexcited dynamics Introducing conduction 
and valence bands within mean field (e.g., effective mass) 
approximation, these interactions can be partitioned into 
two groups P First group combines the interactions which 
conserve the number of electrons/holes, and can be as- 
sociated with the diagonal terms in the electron-hole in- 
teraction Hamiltonian. In other words, these interac- 
tions are acting within exciton bands of different multi- 
plicity, and therefore, will be referred as the intraband 
Coulomb interactions. The intraband interactions deter- 
mine the binding energies of excitons, trions, and multi- 
excitonsP"^ Spectroscopically, these quantities can be 
measured as energy shifts in transient-absorption and flu- 
orescent techniques. One of the important areas where 
such interactions are widely exploited is wav e-fun ction 
engineering in semiconductor heterostructure d 9 ! 10 ! with 
applications to the problem of the lasing in single-exciton 
regime PI 

Second group includes off-diagonal electron-hole inter- 
action terms in the Hamiltonian. These interactions lead 
to the valence-conduction band transitions which change 
the number of electrons and holes conserving the total 
charge. In other words, the off-diagonal Coulomb in- 
teractions are responsible for the photoexcited dynam- 
ics which involves transitions between exciton bands of 
different multiplicity. Therefore, we refer to such in- 
teractions as the interband Coulomb interactions. Im- 
portant examples involving the interband transitions are 
non-radiative Auger recombination and impact ioniza- 
tion processes!^ 

Lately, carrier multiplication (CM) in semiconductor 
NCs has received significant attention ! 13 " 17 ! Also referred 
as the multi-exciton generation, this is a process of more 
than one electron-hole pair generation per absorbed pho- 
ton which has great potential for applications in photo- 
voltaic, photochemical, and energy storage devicesP^l 
A number of models for CM dynamics based on the inter- 



band Coulomb interactions which go beyond simple im- 
pact ionization model have been proposed recentlyP^U 
Direct experimental technique for measuring interband 
Coulomb interactions in NCs should greatly advance our 
understanding of CM by shedding new light on the multi- 
exciton photogeneration dynamics and providing verifica- 
tion of the proposed theoretical models. 

Since CM takes place in the region of high exciton and 
biexciton DOS, key features associated with a variety 
of the inter- and intraband interactions are spectrally 
overlapped and cannot be resolved using conventional 
techniques, e.g., transient absorption and time-resolved 
fluorescence spectroscopies. In this case, the resolution 
can be tremendously enhanced by using coherent ultra- 
fast nonlinear techniques in which various interaction 
pathways are distinguished through the spatial phase- 
matching conditions Further representing the coher- 
ent signal as a 2D Fourier transform polarization with 
respect to various pulse delay times, one can gain direct 
insight into the hidden interactions by looking at the po- 
sition, line-shape, intensity and phase of the associated 
cross-peaks ! 31 " 35 ^ The latter techniques are referred as 2D 
correlation spectroscopies and represent ultrafst optical 
counterpart of 2D NMRP 

Among 2D correlation spectroscopies, the double- 
quantum coherence technique proposed by Mukamel and 
co-workers is considered to be most sensitive to the many- 
body correlations! 3 - 7 -^ This technique has been applied to 
study the intrerband Coulomb correlations leading to 
biexciton formation in organic molecular systems f * 39 * 
and semiconductor quantum welld^^. The application 
of this methodology to probe interband Coulomb inter- 
actions, requires additional theoretical study which is the 
focus of this paper. 

To model nonlinear optical response of coupled 
exciton-biexciton states, we adopt the exciton scattering 
model previously developed by us to provide a unified 
treatment of CM dynamics in semiconductor NCsP^ To 
avoid difficulties associated with large number of reso- 
nances in the region of high-DOS manifold and to un- 
ambiguously identify 2D spectroscopic signatures of the 
interband interactions, we performed numerical calcula- 



2 



tions using reduced level scheme with the parameters for 
PbSe NCsPl 

The paper is organized as follows: In Sec. [TTJ we dis- 
cuss our general theoretical model for the nonlinear opti- 
cal response of coupled exciton-biexciton states and spec- 
troscopic features immediately following form the theory. 



In Sec.|III[ the numerical calculations of double-quantum 
coherence signal for PbSe NCs are discussed in details. 
Finally, we present our conclusions in Sec. |IV| 



II. 



NONLINEAR RESPONSE FROM COUPLED 
EXCITON AND BIEXCITON STATES 



Let us consider an ensemble of NCs, in which the car- 
rier dynamics is restricted to the coupled exciton and 
biexciton manifolds described by the following projected 
Hamiltonian 29 



H = H + V c . 



(1) 



Here, the first term, 



H = ^\x a )hhj x a {x a \ 



E 



xx k )hhj% x {xx k \, (2) 



describes non-interacting exciton, \x a ), and biexciton, 
\xxk), states characterized by the energies hu x and 
hu^ x y respectively. These energies already include corre- 
sponding binding energies associated with the intraband 
Coulomb interactions. The second term 



Vc = ^^\x a )V*'l x {xx k \+h.c., 



(3) 



represents the interband Coulomb interactions, V x ^ x , be- 
tween the states. In this paper, we drop the vacuum, 
\xq) : to biexciton couplings, V£q ,x , considered in the CM 
theoryP^U since their contribution to the nonlinear op- 
tical response is negligible. Explicit representation for 
interaction matrix elements, V^'j^, in terms of the single- 
particle couplings, and related matrix equations defining 
the exciton and biexciton states can be found in Ref [29 . 

Within the exciton scattering model, coupled exci- 
ton and biexciton state propagation is described by the 
Hilbert space Green function 29 



G(t) 



G x {t) G x ' xx 
G xx ' x {t) G xx 



(*) 
(t) 



(4) 



whose diagonal blocks determine the intraband transi- 
tions renormalized by the even-order interband scattering 
events. Since, these propagators do not change the exci- 
ton state multiplicity, we will refer to the associated pro- 
cesses as intraband propagation. The off-diagonal compo- 
nents describe the odd-order interband scattering events 
changing the exciton multiplicity. As a result, associated 
processes will be referred below as the interband propa- 
gation. A closed set of equations for the Green function 
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FIG. 1: Double-sided Feynman diagrams associated with 
Eq. (|8| for the third-order nonlinear optical polarization char- 
acterized by the wave vector k s = ki + k2 — k3. 



based on the scattering matrix formalism is provided in 
Appendix [Aj 

The interaction of coupled exciton and biexciton states 
with optical field is described by the extended Hamilto- 
nian 



H opt = H — p,E(i), 



(5) 



where H is the material part (Eqs. |l])-([3|). The second 
term in Eq. ([5| contains electric field 



E(r,t)=J2£j(t-T j )e a «- 



r—iuJot 



C.C.. 



(6) 



which is a sequence of four linearly polarized non- 
overlapping ultrafast pulses. The first three pulses excite 
the third order nonlinear polarization and the fourth one 
is used to heterodyne the signal!™ Each pulse is char- 
acterized by a wave vector k^ and central frequency coo 
which is set identical for the whole train. We assume im- 
pulsive excitation regime, and approximate the envelope 
functions as £ ) (t — Ti) = £o5(t — Ti), where S(t — Ti) is the 
Dirac delta- function and £q is effective amplitude. The 
latter without lack of generality is set to £q = 1. 

The dipole moment operator projected on the optical 
field in Eq. ([5| has the following block-matrix form 



/i 



ft ft 



(7) 



where the upper (lower) diagonal block describes the in- 
terband vacuum-exciton transitions and the intraband 
exciton (biexciton) transitions. The off-diagonal blocks 
describe interband exciton-biexciton transitions. 

Using the sum-over-states approach along with the 
rotating-wave approximation™, we derived the follow- 
ing expression for the nonlinear polarization propagating 
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along the k s = ki + k 2 — k 3 direction 
P(*3,*2,*i) = (Tr[A-G(t 3 )A"G(t 2 )A + G(ti)A + p] (8) 
- Tr\fi-G\h)jl-G{U + i 2 )A + G(ti)A + p]) 



x e iu;o^3+2iu;o^2+'iwo^i 

containing introduced Hilbert space Green functions, 
G(t) . In Eq. ([8}, the dipole operator (Eq. Q) is par- 
titioned into the sum of two terms, fi = fi + + de- 
scribing optical transitions up and down, respectively.^ 
p = |#o)(#o| is the equilibrium density operator for the 
ensemble of NCs in the vacuum state. 46 

The first (second) term in Eq. ([8| are associated 
with the Liouville space pathway schematically depicted 
by first (second) double-sided Feynman diagram in 
Fig. ([I]). These pathways have been studied before for 
exciton and biexciton manifolds, where the intraband 
dipole tr ansitions and the interband propagations are 
forbiddenPMEoEIIIl However, in NCs these processes 
should be accounted for and, as shown below, lead to 
new features in double-quantum coherence signal, which 
is defined as 

pOO POO 

S(n 3 ,n 2 ) = / dt s dt 2 P(t 3 ,t 2 ,t 1 = 0) (9) 

J JO 

x e i(^ 3 -o;o)^3+*(^2-2u;o)^2 

Note that for further convenience, the origin of (^3,^2) 
2D signal representation is shifted to (ljq, 2cjo), where uo 
is pulse central frequency. 

To calculate double-quantum coherence signal accord- 
ing to Eqs. ([sj) and one needs to find the Green 
function (Eq. Q) by solving numerically set of equa- 
tions provided in Appendix [Aj However, this general ap- 
proach does not give insight into the cross-peak structure. 
Therefore, we apply a different approach using the follow- 
ing representation for the diagonal blocks of the Hilbert 
space Green function^ 

G xx (t) = J2Y,\ xx *) A ki^ xx i\ e ~ iG,tt - ( n ) 

£ kl 

Here, complex frequencies, Cb^ = cj£ — 27^, are the poles of 
G(cj), representing the eigenfrequencies of the Coulomb 
coupled exciton and biexciton states. Hereafter, we will 
refer to these states, as the quasiparticle states. 
KM) = res{G* 6 (^)} (Aff (w 5 ) = res{G%f(^)}) are 
exciton (biexciton) Green function residues determining 
probability amplitudes for the transition between \x a ) 
and \xb) (\xxk) and \xx{)) states associated with the time 
evolution of a quasiparticle state 

Similar expansion for the Green function off-diagonal 
blocks reads 




FIG. 2: Level diagrams representing optical transition and co- 
herence propagation pathways for the components of double- 
quantum coherence signal given by Eqs. (13) and (14). 



G x,xx (t) 



a,k 



(12) 



Solid arrows represent optically induced transitions while the 
dashed arrows show coherence propagations during pulse de- 
lay times. Notations for the dipole moment elements and 
transition amplitudes are the same as used in Table [I] 



where A^'^(^) = res{G^ x (u^)} are the transition 
probability amplitudes between exciton, \x a ), and biexci- 
ton, \xxk), states associated with Adopted represen- 
tation (Eqs. ( 10 )— ([12])) provides connection between the 
spectroscopic resonances determined by the quasiparticle 
frequencies and exciton/biexciton dynamics described by 
the associated transition amplitudes. 
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TABLE I: List of coefficients associated with the double-quantum coherence signal components given in Eq. (14) 



Using Eqs. ([To|-([T2| along with Eqs. one finds 

that double-quantum coherence signal can be represented 
as a sum of eight components 



5(0 3 ,0 2 ) = ^5^(03,02). 
i=i 



(13) 



Each component in Eq. (13) can be attributed to specific 
transition (Liouville space) pathway shown in level dia- 
grams of Fig. [2] Mathematically, signal components can 
be represented in the following compact form 



a(3 qj 



^3 — ^rv 



UJ{3 



(14) 



where gj^q, = tip — cD*. Expressions for the coefficients 
^^•(^a) an d B qj (Cjp) in terms of the dipole matrix ele- 
ments and the transition amplitudes are listed in Table [l| 
According to Fig. |2j the signal components Si and 62 
depend on the transition pathways which comprise only 
the intraband propagation events following the dipole 
transition and no interband propagation. This shows up 
in Table [T] as the dependence of the associated coefficients 
A qj (ti a ) and B qj (tip) on intraband transition amplitudes 
A^ and A xx only. The interband transition amplitude, 
A x,xx , appears in B qj {tip) (A qj (ti a )) for the signal com- 
ponents 63 and 64 (S5 and Sq) and reflects an inter- 
band propagation event in high, i.e, uop ~ 2co>o, (low, i.e., 
uj a ~ oj ) frequency region (Fig. |2|. Finally, the last two 
components 5V and 5s depend on the pathways which in- 
volve both high and low frequency interband propagation 
events (Fig. [2]). 



Using Eq. (14) and Table [T| a general analysis of the 
cross-peak structure can be performed. For this pur- 
pose, we initially assume that the interband Coulomb 
interaction is turned off. In this case A x,xx = 0, and 
thus, all signal components vanish except Si and 82- 
Further taking into account that A x b (uj^) = 5 a bA x a (u x ) 
(Aj£f (ug) = SkiA xx (uj^ x ))^ entering expressions for Si and 
S2, become dependent on the unperturbed exciton (biex- 
citon) frequencies, the following two groups of (^3,^2)- 
cross-peaks can be identified: (uj x ,uj x ) and (uj x ,uj x b ) due 
to Si (Q3, ^2), and (uj xx ,uj x ) and (uj xx , uj x,xx ) clue to 
S2(^3,^2)- Note, that as a result of the intraband exci- 
ton transitions (Fig. the cross-peaks associated with 
Si depend on the exciton frequency only. 

Turning on the interband Coulomb interactions leads 
to the mixing of the exciton and biexciton states and 
splitting of uj x and uj^ x frequencies into quasiparticle fre- 
quencies uj£. The splitting lead to the appearance of new 
cross-peaks in all components of the signal. Further anal- 
ysis of new cross-peaks requires their identification using 
simplified level structure. Before proceeding with this 
analysis, we notice that in the important weak Coulomb 
limit of the exciton scattering model multiple scattering 
events between coupled states reduce to a single (Born) 
scattering event l 2 ^ Expressions for double-quantum co- 
herence signal in this limit are derived in Appendix [Bl 



III. NUMERICAL CALCULATIONS AND 
DISCUSSION 

To identify unique signatures of the interband interac- 
tions in the 2D spectra, we need to reduce the number 
of cross-peaks by adopting a reduced level scheme shown 
in Fig. [3] (b). Comparison of this model with y x ' xx in- 
teraction added and generic multi-levle scheme (Fig. [5| 



5 



{a) Exciton | v 
band \ x l 



Biexciton i, 




(b) Exciton . v 
band \ x / 



l*2> . 






X 












1 












9 i 


MlO 



FIG. 3: Level diagrams for five- level systems reflecting dif- 
ferent excitation regimes in NCs: (a) ftwo ~ 2E g and (b) 
/kJo > 2^. Solid arrows show optically allowed transitions 
and dashed arrows stay for the interband Coulomb interac- 
tions. Accordingly, the signal associated with panel (a) has 
contribution from Si, &, S3, and S4 pathways shown in 
Fig. [2] The signal associated with panel (b) is contributed 
by Si, &, S5, and Sq pathways shown in Fig. [2] 



shows that the former one captures all the transition 
pathways present in the general case. As a result, double- 
quantum coherence signal from the five-level system also 
has eight non- vanishing components explic itly c alcul ated 
in Appendix [C| Obtained expressions (Eq. (C4)-(C6) and 
Table [ll| are used in our numerical calculations. The 
analysis of the signal can be further simplified by con- 
sidering two excitation regimes, schematically shown in 
panels (a) and (b) of Fig. [3] These regimes are based on 
specific properties of exciton and biexciton level struc- 
ture in semiconductor NCs. They also allow one to elim- 
inate the contributions of S7 and Ss which include two 
interband propagation events (Fig. [2| complicating the 
cross-peak patterns. 



We assume that considered semiconductor NCs pos- 
sess inversion symmetry leading to the parity conserva- 
tion. We further assume that the lowest exciton state 
is optically active ("bright") and, thus, has odd par- 
ity. Since, the lowest biexciton state is formed from the 
two lowest exciton states, it has even parity. The inter- 
band Coulomb interaction conserves parity and, there- 
fore, does not couple exciton and biexciton states with 
different parities. Such an exciton and biexciton level 
structure is realized in lead-salt semiconductors. There- 
fore, without lack of generality, we adopt the effective 
mass model based on the k • p-Hamiltonian for PbSe 
NCsP The diameter of the NCs is set to 5 nm. Since 
the calculations of the exciton and biexciton binding en- 
ergies, as well as the interband Coulomb matrix elements 
are not directly accessible from the adopted model, we 
will consider these quantities as parameters. To explore 
2D signatures of the interband interactions, the interac- 
tions strengths will be varied in a broad range. 



A. Excitation Regime ftwo ~ E g . 

First, we consider the case shown in panel (a) of Fig. 3j 
in which pulse central frequency is resonant with the 
lowest exciton state, \x\), i.e., Huoq ~ E g . Under this 
condition, the double excitation energy corresponds to 
the Coulomb coupled even-parity lowest biexciton state, 
I ££2) and even-parity exciton state \x2). These sates are 
allowed for optical transitions from \xi). According to 
the adopted effective mass model, we set Hoof = 1.0 eV, 
hu x = 2.2 eV, and hu xx = 1.8 eV with the bind- 
ing energies included. Noteworthy, this model is sim- 
ilar to the so-called coherent superposition model in- 
troduced by Shabaev, Efros, and Nozik to study CM 
in PbSe nanocrystals. 27 Since, the interband transition 
dipoles fif and \l^ x depend on the same valence-to- 
conduction band matrix elements, we set them equal. 
The intraband exciton transition dipole, /ifi, vanishes 
in bulk semiconductors due to the quasimomentum con- 
servation constraint. In NCs, this constraint is relaxed 
and the intraband transitions can acquire some oscillator 
strength. Therefore, we consider three different cases in 

which /ifi = 0, /ifi = 0.5/ifo? an d M21 = Mio- I n g enera ^ 
double-quantum coherence signal associated with so de- 
fined level scheme is a superposition of Si, Si, S3, and 
S4 components (Appendix |C| and Fig. [2|. 

To start, we set /ifi = reducing the total signal to 
the S2 component which is associated with the pathway 
involving no interband propagation (Fig. I2J. The corre- 
sponding 2D spectra are shown in column (a) of Fig. Q . 
For y x ^ xx = (upper panel), the 2D spectrum contains 
two (O3, ^2)-cross-peaks. Namely, peak 1 at (u x ,u xx ) 
and peak 2 at (uj xx { x ,uj xx ). The difference between the 
^3-resonances in this case provides the biexciton bind- 
ing energy defined as u xx — 2cjf. The middle and lower 
panels in column (a) show the same spectra for finite in- 
terband Coulomb coupling strength y x ' xx = 0.2 eV and 
yx,xx _ q 5 e y^ respectively. Switching on the inter- 
action leads to the renormalization of the peak 1 and 
2 frequencies as (uf,U2-) and (0J2- —^1,^2-), respec- 
tively. In addition cross-peaks V at (cjf ,6^2+) and 2' at 
(c^2+ — , c^2+) emerge. Here, co>2± denotes mixed (quasi- 
part icle) energy of coupled \xi) and \xx2) states given by 
Eq. (C5). Accordingly, the increase in the interband cou- 



pling leads to the increase in the energy splitting between 
the cross-peaks as observed in Fig. [4j 

Setting /ifi to finite values at V 2 XjXX = (upper panels 
in column (b) and (c) of Fig. [4|, adds Si contribution 
to the 2D signal which is solely associated with the in- 
traband exciton transitions (Fig. |2j. This leads to the 
appearance of already seen cross-peaks V and 2' ', whose 
intensities increase as /ifi increases. Adding the inter- 
band Coulomb interaction (middle and lower panels in 
column (b) and (c)) adds the contributions from the S3 
and S4 signal components associated with the interband 
propagation (Fig. However, the resonances associ- 
ated with these components overlap with already existing 
cross-peaks 1, 2, and V 2! . 
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FIG. 4: The absolute value of double-quantum coherence signal for the excitation regime huuo « in NCs shown schematically 
in Fig. [3] (a). The interband exciton dipole is set to (a) /ifi = 0, (b) /ifi = 0.5/ii , and (c) /ifi = Mio- Within each column 
the interband Coulomb interaction increases from top to bottom taking the following values: y x > xx — Q.0, y x,xx — Q.2 eV and 



V* 



0.5 eV. Red line stands for ^2 = 2^3 



The interference of all pathways contributing to the 
cross-peaks can lead to dramatic variation in their inten- 
sities. As observed in middle and lower panels of col- 
umn (b), the increase of the Coulomb coupling reduces 
the intensities of cross-peaks 1 and 2. For equal strength 
of inter- and intraband transition dipoles (column (c)) 
and strong interband coupling (lower panel), cross-peaks 
1, 2 vanish. This effect can be rationalized through the 
formation of symmetric, |2+) = l/y/2(\x2) + |##2)) an d 
antisymmetric, \2— ) = l/y/2(\x2) — \xx2)), superposi- 
tions of the coupled states where the latter one is opti- 
cally "dark". 

Finally, we note that in the excitation regime Hujq « 
Eg, the contributions of signal components S\ through 
64 correspond to the cross-peak patterns that do not al- 
low one to qualitatively distinguish between the effect of 
the interband Coulomb interactions associated with the 
components 63 and #4 (Fig. [2| and the intraband dipole 
transitions accounted for in Si (Fig. Q). This conclu- 
sion immediately follows from the comparison of the mid- 
dle panel in column (a) and the upper one in column (b), 
which have the same cross-peak patterns. Therefore, con- 
sidered excitation regime does not allow one to obtain a 



clear signatures of the interband Coulomb interactions. 
As we show next, the situation is quite different for the 
excitation regime shown in Fig. [3] (b). 



B. Excitation Regime Huoq > 2E g . 

According to the adopted effective mass model, the 
lowest odd-parity biexciton state in NCs is second-to- 
lowest biexciton state whose energy is slightly higher than 
2E g . This biexciton is Coulomb coupled to an exciton 
state of the same parity. Since the latter exciton state is 
optically allowed, we turn pulse central frequency in reso- 
nance with this state. This leads to the excitation regime 
shown in Fig. [3] (b). Following, the notations used here 
for the five-level system, we denote the coupled exciton 
and biexciton states as \xi) and \xxi), respectively. Solv- 
ing the effective mass model for PbSe NCs^ as well as 
adding the binding energies, we set hoof = 2.2 eV and 
hui x = 2.5 eV. We further chose a pair of double-excited 
even-parity exciton, \x2) and biexciton, \xx2), states with 
the energies fkjf = 3.6 eV and fruf^ =4.1 eV, respec- 
tively. These states can be coupled through the interband 
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FIG. 5: The absolute value of double-quantum coherence 
signal for the excitation regime huo > 2E g in NCs shown 
schematically in Fig. [31(b). The interband transition dipoles 
are set to (a) = M21 = 0, (b) ^1 = M21 = Mio- Within 
each column the interband Coulomb interaction increases 
from top to bottom taking the following values: y x,xx = 0.0, 
y x, xx = o.l eV and V^ xx = 0.4 eV. Red line stands for 

Q 2 = 2Q 3 - 



along Q 2 , since y x ' xx = 0. Instead, the interband inter- 
action, Vi' xx , renormalizes 1,2 cross-peak frequencies as 
(uji-,uj xx ) and (uj xx — uji-,uj xx ) and leads to new fea- 
tures 3 and 4 emerging at the frequencies (oji+,oj xx ) and 
(u xx — cji+, uo xx ), respectively. 

Switching on the intraband optical transitions adds a 
new set of resonances at ^2 = w x as can ^ e seen m 
Fig. [5] (b). Specifically, the cross-peaks 1', 2', 3', and 
4' emerge at {001-, cjf), (w x — uoi-,uo x ), (^1+^2) an d 
(cjf — ^1+7 w x )? respectively. These resonances result from 
the contributions of Si, S5, and $6 signal components. 
These components also provide additional contributions 
to cross-peaks 1 — 4. As the Coulomb interaction rises 
from top to bottom of column (b), the interference effect 
similar to that in Fig. [4](b) and (c), affects the intensity 
of the cross-peaks. 

Comparison of columns (a) and (b) in Fig. Q, shows 
that the adopted excitation regime allows one to distin- 
guish between the effects of the interband Coulomb in- 
teraction and intraband optical transitions: If y x > xx = 
and /ifi = M21 = 0? one observes a pair of cross-peaks 
1, 2 along ^3. The peaks positions are symmetric with 
respect to the line 2^3 — ^2=0 (red in Fig. [4]and[5|, 
and the splitting provides the biexciton binding energy. 
For Vi' xx = and /ifi = M21 ^ 0> a new symmetric 
pair of cross-peaks V and 2' shows up at Q2 = ^f- For 
yxx _^ q anc j jjx^ _ ^xx _ |- wo p a i rs Q f cross-peaks 

(i.e. 1, 2 and 3, 4) symmetric with respect to the line 
2^3 — ^2=0 are observed. In this case the energy split- 
ting between symmetric peaks is due to the biexciton 
binding energy renormalized by the interband Coulomb 
coupling. Finally, for V xx 7^ and /ifi = M21 7^ a ^ 
eight cross-peaks (i.e. 1 — 4 and V — 4') show up in the 
2D spectrum. Note, that the symmetry of the cross-peak 
pairs with respect to 2^3 — ^2 = line can be used to de- 
termine the biexciton binding energies and the interband 
Coulomb interactions in more realistic situation, when 
large number of the resonances make their identification 
difficult. 



Coulomb interactions. However, they belong to the high- 
energy region where these interactions are much weaker 
than y x ' xx and therefore we set y x ' xx = 0. Similar to 
Sec. ~ 



the interband transition dipoles /if , \i 2 



III A, 

and \i x ' 2 are set identical, and the intraband ones are 

varied as /if 2 = /iff = ano ^ M12 = Ml 2 = Mio- ^ n § en ~ 
eral, double-quantum coherence signal associated with 
such defined level scheme (Fig. [3] a) is a superposition of 
Si, S2, S§, and Sq components (Appendix |C| and Fig. [2]). 

Fig. [5] (a) shows double-quantum coherence spectra in 
the case when the intraband optical transitions are for- 
bidden, i.e., /ifi = M21 = 0- ^ interband Coulomb inter- 
action is set to y x ' xx = (upper panel in column (a)), 
the spectrum is solely due to 62 component and con- 



tains cross-peaks 1 and 2 at (ujf,uj xx ) and {uj 2 { 



xx,x %X X 



), 



respectively. This is the same situation as we discussed 
in Sec. |III A| In contrast to Sec. Ill A[ turning on the 
interband interaction does not lead to the level splitting 



IV. CONCLUSIONS 

Using previously developed exciton scattering models 
accounting for the Coulomb interactions between exciton 
and biexciton states, we have derived closed expressions 
for 2D double-quantum coherence signal. In general, the 
Liouville space pathways contributing to the signal ac- 
count for the interband scattering processes and the ef- 
fect of the inter- and intraband optical transitions. The 
interplay of these two effects is important to take into 
account considering spectroscopic probe of photoexcited 
carrier dynamics in semiconductor NCs. Presence of the 
intraband optical transitions and the interband scatter- 
ing effect corresponds to the formation of new cross-peaks 
in the 2D spectra. To find clear spectroscopic signatures 
of each of the effects, we have considered two excitation 
regimes and found that the last one, i.e. Hljq > 2E g 
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carries desired sensitivity. 

Although, the calculations performed here used effec- 
tive mass model for PbSe NCs they can be considered 
rather model since unrealistic values for the biexciton 
binding energies and interband Coulomb couplings have 
been used. The analysis of realistic spectra would re- 
quire more computational efforts and necessity to deal 
with large number of cross-peaks. However, the trends 
we have identified in the 2D spectra using our model 
system should be present and eventually sough in more 
realistic calculations and in the interpretation of experi- 
mental observations. 
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Appendix A: The Exciton Scattering Formalism 



For the projected Hamiltonian given by Eqs. 
the matrix elements of exciton (n = x) and biexciton 
(n = xx) free propagator have the following form 

OJ~UJ k 

where the complex frequency, = — contains 
the fc-th frequency, cjJ?, from the projected Hamiltonian 
(Eq. Q), and the dephasing rate, 7^. The exciton (n = 
x) and biexciton (n = xx) scattering matrix elements 
satisfy a set of linear equations 

E Mi " i<4Mff&M] W = i4M, (A7) 

kl 

with the self-energy matrix elements defined as 

4M = ih- 2 E v:fgU")VT, (A8) 

kl 

where fa = xx (m = x) if n = x (n = xx). Finally, 
the interband scattering matrix elements can be obtained 
from the following linear transformation: 

mn 



In this Appendix, we provide a closed set of Eqs. ( Al ) 



(A9) allowing one to calculate the Hilbert space Green 
function (Eq. Q) using the scattering matrix formalism. 
Derivation details can be found in Ref. [29] . 

We start with the Fourier representation of the time- 
domain Green function 



G(t) 



— G(uj) exp (—loot). 



in terms of the frequency domain counterpart 

1 ; \G xx > x {lo) G xx (lu) J' 
This quantity satisfies the following equation 

G(Lj)=g(u>)+g(Lj)f(u>)g(u>), 
where g(uu) is the intraband free propagator 



g x (u) 

g xx (u) 



and T(uj) is the scattering operator 

f (u) = ( T x {uj) T x ' xx (oj) 
K ' \f xx > x (uj) T xx (uj) 



(Af) 



(A2) 



(A3) 



(A4) 



(A5) 



containing exciton scattering, T x , biexciton scattering, 
T xx , and the interband scattering, T x ' xx , components. 



Appendix B: 2D Signal in Weak Coulomb Limit 

In this Appendix, we use weak Coulomb limit of the 
exciton scattering model 29 to derive the double-quantum 
coherence signal. Adopted limiting case assumes that 
the interband Coulomb coupling is weak compared to 
the energy difference between coupled states and/or their 
broadening. This ultimately leads to the leading contri- 
bution of the interband Born (single) scattering event. 

Ret ainin g only leading Born contributions in 
Eqs. ( |A1| )-(|A9D, one finds that the diagonal blocks 
in Eq. (El) ar<P^ 



G x {t) = Y,W){x a \e-^\ (Bf) 

a>l 

G xx (t) = ^|xx,)(xx,|e-^ £ . (B2) 



k>l 



xx 
k 



Here, the complex frequencies Co x = u x — ij x and Co 
00 k X ~ ^k X con tain exciton, 7^, and biexciton, 7^, de- 
phasing rates, respectively. 

In the adopted limit, the leading off-diagonal blocks in 
Eq. (|4| become^ 1 



&' xx ® = EEw A sr< 



XX k 



(B3) 



x (e 



a>l k>l 



-iCot x t 



rt ), 



where the interband transition amplitude is 



A : 



a.k 



(B4) 
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According to Eqs. Q-([9|, and the Green function rep- 
resentation given by Eqs. (|10|)-( 12 ), double-quantum co- 
herence signal becomes a sum of six components 



s(n 3 ,n 2 ) = X) ^(03,02), 



(B5) 



where each component is a limiting value of the general 
signal components given by Eq. ( 14 ) and can be identified 
with the corresponding pathways in Fig. [2] Terms S7 and 
$8 are neglected, since they have second order in V x,xx 
contributions which are beyond the Born approximation. 

The first two signal components provide zero-order in 
yx,xx contributions 



and 

S[ 1] (fi 3 ,fi 2 ) 



ac,nm 
UJ X)XX 



X \X,XX XX XX, x x /-pi o\ 
A^Oc^ c, m^mn^na l^aO V- D1Z / 



+ (fia-^Kfia-o;-))' 

These terms also provide first order corrections but con- 
tain new cross-peaks associated with the interband scat- 
tering. Specifically, these new peaks appear at ^3 = cD^f , 

^3 =^6,'m' and ^3 =^nm« 



Sl(ft3,ft 2 ) = E 



ft* 



u 6c 



with cD^ c = 



cD^*, and 



5 2 (^3,^2) = E 



^2 - 



n 3 



with u xx,x — u xx 



Next two terms, i.e., Sj with j = 3,4 



S,-(fi 3 ,fi2) 



E 

bc,n 



F 3 

bc,n 



n 2 



1 



^3 - 



ft. 



bc,n 



UJ 



ft 3 



ft 3 



~ XX, X 

W n,c 



with 



abc,n 
abc,n 



X X,XX A XX,X X X 

H'OcH'cn iV n, b^ba^aO 

x x \x,xx xx,x x 
^Ocrcb^b, n V n, oPaO' 



(B6) 



(B7) 



(B8) 



(B9) 
(BIO) 



provide first order corrections to Eq. (B6) and Eq. (|B7 



Note that the resonances in Eq. (B8) coincide with the 
resonances in Eq. (B6) and Eq. (B7). 



The last two contributions to the signal read 

s 3 (n 3 ,n 2 ) = I £ ^ c a x c ; x X'>L^ (bh) 

1 



abc,m 
, .x,xx 



ft 2 



n 3 
1 



•W?)(ft3-W*f) 



(ft 3 -^ c )(ft 3 



Appendix C: 2D Signal associated with Five-Level 
System 

In this Appendix, we derive the expressions for double- 
quantum coherence signal associated with five-level sys- 
tem shown in Fig. [3] (b). For completeness, we also add 
the interband Coulomb interaction y^ xx i n the high- 
energy region. 

In this case, the non-interacting exciton and biexciton 
components of the projected Hamiltonian (Eqs. (|2|-([3|) 
reduce to the following form 



#0 



E 

a=l,2 



x a )hh:l{x a \ + E \xx k )hul x (xx k \.(Cl) 



k=l,2 



Here, the exciton energies are set to be hoof ~ ujo and 
hw x ~ 2uoq. The difference between adjacent exciton and 
biexciton energies is assumed to be within the pulse spec- 
tral widths a <C c^o, i.e., \uo x — uo xx \ < a for p = 1, 2. The 
adjacent exciton and biexciton states are coupled by the 
interband Coulomb interaction 

Vc = E \xc)V?> xx (xx c \ + h.c. (C2) 

c=l,2 

Finally, the transition dipole operator for the adopted 
five-level system reads 



A = ko)w5i(^i| + \x 1 )^ x 2 (x 2 \ 

+ kl)Ml,'T(^2| + |xXi)/i^'2 (X 2 \ 
+ \xXi)fl X 2 (x%2\ + ^- c -- 



(C3) 



For so defined model, all eight double-quantum coher- 
ence signal components (Eq. (14) and Table [i]), reduce to 
the following form 



a,(3=± 



o 2 



(C4) 



4* 
^3 



where oj 2 /3,ia 
are 



: ^ 2 /3 ' 



Col a1 and the quasiparticle energies 
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Number of 
Interband 
Propa- 
gation 
Events 


Si 




Aj 







Si 

s 2 


Im§iI 2 ImT 2 | 2 
\»oi\ 2 \» x G\ 2 


Af« 
Af„ 


A XX 

A 2/3 




s 3 


I^SilV?;r^i 


Af„ 


A XX, X 

A 2/3 


1 


s 4 


Im§iI 2 mIiM2;T 


Ala 


A 2£ 




s 5 


i x |2 xx,x x 
iMoil A* 1, 2^21 


A x,xx 


AX 

A 2/3 




S 6 


iMoi | 2 A*f,2A*"'r 


A. 


A XX 

^*-2/3 


2 


S 7 

s 8 


KiH/^l 2 
Ia&IVi™! 2 


A x,xx 
A la 

a x,xx 
A lc* 


A XX, X 

A 2/3 

A x,xx 
A 2/3 



TABLE II: List of coefficients associated with the double- 



quantum coherence signal components given in Eq. (C4) 



with c = 1,2. The expressions for the coefficients Aj 
and Bj are listed in Table [IT] and contain the following 
transition amplitudes^ 



A XX 
A C± 



AX, XX 
c± 



(uj c± - u xx ) 



(^c+ - 




)' 


(^c± - 




) 






)' 


K c 







h(u>, 



c+ 



-)' 



(C6) 
(C7) 
(C8) 



with c = 1, 2. 



Sec. Ill A 



To obtain the signal considered in 
(Fig. [3] (a)), one has to set y x ' xx = 0. This reduces 
contributions to the sum of Si, S2, S3, and S4 compo- 
nents. The signal considered in Sec. Ill B| (Fig. [3] (b)) 
can be obtained by setting y x ' xx = 0. This results in the 
following non- vanishing contributions: Si, S2, S5, and 
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